Theoretical studies of spin-dependent electronic transport in ferromagnetically 

contacted graphene flakes 
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Based on a tight-binding model and a recursive Green's function technique, spin-depentent bal- 
listic transport through tinny graphene sheets (flakes) is studied. The main interest is focussed 
on: electrical conductivity, giant magnetoresistance (GMR) and shot noise. It is shown that when 
graphene flakes are sandwiched between two ferromagnetic electrodes, the resulting GMR coeffi- 
cient may be quite significant. This statement holds true both for zigzag and armchair chiralities, 
as well as for different aspect (width/length) ratios. Remarkably, in absolute values the GMR of the 
armchair-edge graphene flakes is systematically greater than that corresponding to the zigzag-edge 
graphene flakes. This finding is attributed to the different degree of conduction channel mixing 
for the two chiralities in question. It is also shown that for big aspect ratio flakes, 3-dimensional 
end-contacted leads, very much like invasive contacts, result in non-universal behavior of both con- 
ductivity and Fano factor. 

PACS numbers: 81.05.Uw, 75.47.De, 75.47.Jn 



I. INTRODUCTION 



Recently a lot of interest has been directed to carbon- 
based systems, in search for alternative materials which 
would make it possible to go beyond the silicon tech- 
nology. While looking at history of studies along this 
line, one can point out two important milestones: (i) dis- 
covery of carbon nanotubes (CNTs), honeycomb-lattice 
cylinders, dating back to 1991^ and (ii) fabrication of 
individual atomic planes, called graphene (Gr), by exfoli- 
ation from graphite in 2004^ So far, on obvious grounds, 
the CNTs have been much more thoroughly investigated 
than Gr, but this difference diminishes very quickly. 
This paper focusses on spin transport problems, mostly 
on giant magnetoresistance, related to potential appli- 
cations of graphitic nanostructurcs in spintronics. In 
this respect, quite a lot has been done in the case of 
CNTs, there are hundreds of both experimental^^i 
and theoretical ^ 10 ' 11 papers, covering all the transport 
regimes (ballistic, Coulomb blockade, Schottky barrier 
and Kondo) which show that the GMR or TMR (T for 
tunnel) effects are usually quite considerable. The re- 
spective studies on Gr are still scarce. The pioneering ex- 
perimental paper on magnetoresistance is Ref(l2l] , where 
Gr spin valve devices with permalloy contacts have been 
shown to have the GMR effect ~ 10% at room tem- 
perature. The following experimental paper,— reports 
on conductance of Gr showing Fabry-Perot-like patterns, 
and pronounced oscillations of GMR (including changes 
in sign). From the theoretical point of view, the debate 
is still on, and the question whether or not the Gr-based 
spin valves have a good performance is still open. The 
results published so far show that the answer to this ques- 
tion depends critically on the contacts (cf. Ref.flij and 
Refs. jl5lJl6| ). In this study, in contrast to those reported 



hitherto by other theoreticians, a 3-dimensional contacts 
are used. 

The paper is organized as follows: In Sec. II the the- 
oretical method based on a tight-binding model, as well 
as the way the graphene sheets and the contacts have 
been modelled, are shortly outlined. Sec. Ill is devoted 
to presentation of results, whereas the subsequent section 
summarizes the main results. 



II. METHOD AND MODELLING 



The method employed here is similar to that used ear- 
lier while dealing with carbon nanotubes inserted be- 
tween ferromagnetic contacts ^ 10 ' 11 except that carbon 
allotrope of interest now is graphene . The single tt- 
orbital tight-binding Hamiltonian describes a graphene 
sheet of width W and length L (in the current direc- 
tion). The semi-infinite metallic electrodes extend from 
< —oo,0 > and < L,oo > for the source and drain, re- 
spectively. The ferromagnetic electrodes are supposed to 
have spin-split s-bands, mimicking d-bands of real tran- 
sition metals. The total Hamiltonian reads 



H 



ti,j \i,(x}(a,j\ + 



i,a (a, i 



(1) 



where i and j run over the whole device (i.e graphene 
and the electrodes), a is the spin index, and tij and ej i(T 
stand for the hopping integrals and the on-site potentials, 
respectively. The systems under study are impurity-free 
with well transparent interfaces (strong coupling limit) so 
neglecting of correlations in the Hamiltonian should be 
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functions in k-space. 
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FIG. 1: (Color online) Graphene sheets attached to semi- 
infinite 3-dimensional metallic electrodes. The left and right 
parts correspond to the zigzag- and armchair-edge ribbons, 
respectively. The lower part shows details of the interfaces, 
whereas the letters A, B denote a type of the sublattice the 
interface carbon atoms belong to. 
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where a — \f2acr is the fcc-metal lattice constant, g a is 
the surface Green function for spin a, and A a is a rigid 
band splitting chosen so as to give a desirable spin po- 
larization P of the electrodes. Here P=50% has been 



set, corresponding to A 



-2.32 1 and A 



1.6* 



(in the paramagnetic case A CT = — 0.86 i). Inciden- 
tally, this simple paramctrization has already been shown 
to work satisfactorily well in the case of carbon nan- 
otube/ferromagnet systemsi2»i£ £er(k) in Eq. (2) is the 
fcc(lll)-surface energy spectrum, calculated for a semi- 
infinite metal slab according to the method described in 
Ref. The surface k- vectors lie on the metal surface 
but upon attachment of graphene they arc no longer good 
quantum numbers, so the trick is to perform a Fourier 
transformation to the real space and work with those 
surface Green function which are close to the graphene 
interface. After having transformed the surface Green 
function to the real-space, the self-energies and the 
corresponding spectral functions arc computed from 
= Tg%T^ and = i(££ - E^ f ), respectively. With 
a =L or R, referring the left and right electrodes, and 
T being the Gr/electrode coupling matrices. Henceforth, 
the spin indexes a will be skipped for brevity. 



The recursive method goes as follow; 
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justified if one restricts oneself just to the ballistic trans- 
port regime, putting aside the Coulomb blockade^^ and 
Kondo^ regimes. In this regard, the on-site parameters 
in the present approach are used to take into account the 
effect of the gate voltage in the graphene sheet as well as 
the spin band-splitting in the metallic electrodes. 

Fig. Q]shows schematically devices of the present inter- 
est, viz. zigzag-edge graphene (zz-Gr) and armchair-edge 
graphene (ac-Gr). The graphene sheets (black spheres) 
are sandwiched between metallic contacts (light spheres). 
The diameters of the spheres correspond to the nearest 
neighbor spacings. The unit cells in the vertical direction 
are the blunt saw-teeth lines for zz-Gr, and double zigzag 
lines for ac-Gr. As readily seen, it is assumed that there 
is a perfect lattice matching between graphene and the 
electrodes. In fact this assumption is acceptable since the 
graphene lattice constant (ac r = 2.46A) fits really well 
to interatomic distances in such metals like: Cu (2.51 A), 
Ni (2.49 A) or Co (2.55 A)J£ 

The use of a single-orbital electrodes allows to write 
an analytic expression for the electrodes surface Green 



<? L (0) = g L , g B ,{N + 1) = g Rl 
9l,r{%) = {E-Di-^L^i))- 1 

E_r(i) = T iti+ ig R (i + 1)T 4+M , (3) 

G i = {E-D i -E L (i)-E R (i))- 1 . (4) 

Above, gL,R are local Green's functions for the i-th unit 
cell of graphene, the matrices D and T stand for the diag- 
onal and off-diagonal Hamiltonian sub-matrices, whereas 
the full Green's function is given by Eq.fJ] The graphene 
unit cells extends from i=l (following the L electrode) 
up to i=N (preceding the R electrode). So the recursion 
starts with the metal- interface Green functions c/l(0) and 
9r[N+1) being Fourier components of the surface Green 
functions defined in Eq.®. 

The other quantities of interest are transmission (T), 
conductivity ((7) , and shot noise Fano factor (F), as 
well as giant magnetoresistance (GMR). In the ballistic 
transport regime and at zero temperature these quanti- 
ties read: 
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FIG. 2: (Color online) Giant magnetoresistance for ac-Gr 
(solid curve) and zz-Gr (dash curve) flakes vs energy, and 
for 3 different aspect ratios. Note that the number of quasi- 
periods is roughly inverse proportional to A (i.e. proportional 
to the length), meaning that Fabry-Perot-like resonances take 
place. 



FIG. 3: (Color online) Shot noise Fano factor for the armchair 
graphene ca 8nm in width and with the aspect ratio A=0.36: 
for ferromagnetic electrodes (parallel alignment - thin solid 
line, and antiparallel alignment - dash line). Inset presents 
the respective conductivities. 



1 — I i Lril i <_Tj , 

a = (L/W)jTr[T(E F )], 

F = Tr[T{E){l-T{E))]/Tr[T{E)l 
GMR = 100(1 -<7 Ti4 ./(7 T>T ), (5) 

where the arrows tt and tl denote parallel and an- 
tiparallel alignments of ferromagnetic electrodes. 



III. RESULTS 

We use Eqsj5]to determine the quantities of our main 
interest. At first, for comparison, we have calculated 
GMR for zz- and ac-Gr sheets vs. energy which, in princi- 
ple, is proportional to the gate voltage. From the results 
in Fig[3]one can see that GMR, at least in the ac-Gr case, 
is quasi-periodic with the period-length roughly propor- 
tional to the aspect ratio A, i.c inverse proportional to the 
length of the graphene sheets L. Numerically, the periods 
are close to the well-established value p = 2eV /L[nm], 21 
or in the present units p = 0.8 A/W[nm] (the energy unit 
is t=2.7 cV). Incidentally, some features characteristic for 
Fabry-Pcrot resonances in graphene have been already 
reporte d 21 i 22 , where it has been also stressed that irregu- 
larities and defects of graphene-edges may turn the con- 



ventional Fabry-Perot picture into a spectacular quan- 
tum billiard-type one. It is seen in Fig. [5] that the GMR 
factors for ac-Gr and zz-Gr differ in 3 respects: (i) mag- 
nitude, (ii) the maxima of the former are roughly equidis- 
tant in contrast to the latter, moreover (iii) in the former 
the onset of the second sub-band is visible (at E = 0.08, 
except for A=3.4, where W > L), whereas in the latter it 
always is washed out. Qualitatively these results, show- 
ing that GMR vs. energy (gate voltage) changes also in 
sign, are consistent with the experiment of Ref. [l3T ] . 

Another interesting point encountered here, is the uni- 
versal behavior issue due to evanescent modes at the 
Dirac points. It is well-known that graphene sheets of 
large sizes, with W — > oo and L — > oo (with W > L), 
if homogenous, reveal universal values of (7 = 4/ir e 2 /h 
and F = 1/3^^-. The term "homogenous" in this 
context means that the device at hand is an all-carbon 
system, typically, with a central graphene sheet, and 
highly doped graphene electrodes. This theoretical con- 
cept might be quite realistic for widely applied experi- 
mental setups with side- contacted grahene sheets, pro- 
vided evaporated metal contacts are not invasive, i.e. 
they do not destroy the underlaying honey-comb lat- 
tice but merely dope it slightly. It has been also shown 
theoretically that in the case of end- contacted devices 
to 2-dimensional square-lattice contacts, the situation is 
more complicated and respective values of the conduc- 
tivity and the Fano factor do depend on the on-site po- 
tentials in electrodes and on details of bonds between 
square-lattice electrodes and graphcnc<21 Importance of 
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the interface conditions which determine, inter alia, mix- 
ing of propagating modes, and whether or not the mo- 
mentum component along the interface direction is con- 
served, has been noticed in Ref . [l^mmiOIIllIil ' lt 
has been also demonstrated using one-parameter scaling 
arguments that at the Dirac points, conductivity of an 
infinitely large graphene with disorder is infinite (zero) 
if there is not (there is) intervalley scattering^ Beyond 
the Dirac points the evanescent modes give way to the 
propagating ones, and in the ballistic transport regime, 
perfect noiseless transmission due to Fabry-Perot reso- 
nances may take placei 13 i 24 

Fig[3] clearly shows that this scenario also holds in the 
ferromagnetic case for ac-Gr if the aspect ratio is not 
too big (Gr ribbon is long enough). Indeed maxima in 
(7 (Inset) correspond to minima of F (at F=0) in the 
main panel. In the zz-Gr case with the same A value the 
situation is similar (not plotted), but then the oscillations 
in the P and AP configurations are not phase-shifted. 

As regards Gr sheets with bigger A, FigsfJ] and [S] show 
some features of F and G which happen to be similar 
to the universal ones when the Gr sheets are paramag- 
netically contacted (thick black lines). However, if the 
contacts are ferromagnetic there is a tendency for F to in- 
crease and for O - to decrease. This suggests that in gen- 
eral Gr-flakes studied here show non-universal behavior, 
which is to be attributed to the finite sizes, and the inho- 
mogencity resulting in contact-dependent charge trans- 
fer between interface atoms (and accompanying electron- 
hole asymmetry). Incidentally, the present approach 
was shown earlier to yield short-range charge transfer, 
affecting mainly interface carbon monolayers! 10 ' 20 As a 
matter of fact, the non-universal behavior of F and <7 
in Gr-systems has been already reported to be due to: 
charged impurity scattering,— strong disorder—, invasive 
contactsi2£>22 

Finally, it should be noted that the results presented 
in Figs. 2-5 hardly depend on the hopping parameters 
as long as the interfaces are transparent enough. In the 
present theory this condition is fulfilled provided the hop- 
ping parameter across the interface (t c ) is not too differ- 
ent from the geometric mean of hopping parameters for 
the metal electrodes (£m) and graphene (t), i.e. from 
t c = \/tUt- Otherwise, in case of a drastic differentiation 
of the hopping parameters, and opaque interfaces (very 
small t c ) the Coulomb blockade physics may come into 
play?i£ out of reach of the present approach. 



IV. CONCLUSION 

Summarizing, the aim of this study has been to esti- 
mate the effect of the chirality as well as the aspect ratio 
on the GMR coefficient of graphene flakes end-contacted 
to ferromagnetic electrodes. In contrast to other theo- 
retical approaches, the leads are not supposed here to be 
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FIG. 4: (Color online) Shot noise Fano factor for the zigzag- 
edge graphene ca 8nm in width and with the aspect ratio 
A=3.4: for paramagnetic electrodes (thick solid line), ferro- 
magnetic electrodes (parallel alignment - thin solid line, and 
antiparallel alignment - dash line). Inset presents the respec- 
tive conductivities. 
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FIG. 5: (Color online) As|4]but for armchair-edge graphene. 



two-dimensional (of honeycomb- or square-lattice type), 
but more realistically they are modelled as 3-dimensional 
fcc-(lll) semi-infinite slabs. It turns out that for long 
and narrow systems (small aspect ratio), the GMR is 
a quasi-periodic function of energy (gate voltage) with 
the period roughly proportional to the aspect ratio, re- 
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fleeting Fabry-Pcrot-likc resonances, typical of ballistic 
transport. Notably, the GMR coefficient of the ac-Gr 
flakes may exceed 20% - 40%, depending on the aspect 
ratio value, whereas for zz-Gr flakes the corresponding 
figures are distinctly smaller, but still significant. The 
difference is due to the fact that in the case of armchair- 
edge sheets all the interface carbon atoms belong to the 
same sub-lattice (say, A- type), as opposed to the zigzag 
ones with interface carbon atoms of both A- and B-typc. 
This inevitably facilitates intervalley mode mixing in the 
latter case. Another noteworthy point, in the context 
of zero-energy (Dirac-point) conductivity and Fano fac- 
tor, is that for the big aspect ratio. Such systems stud- 
ied here show non-universal behavior when ferromagnetic 
electrodes are applied, ft is so even if the corresponding 
values for the paramagnetic electrodes happen to be close 
to the universal ones (4/7T e 2 /h and f/3 for (7 and F, re- 
spectively). However, these systems are inhomogeneous 
and rather small, so neither the limits W — > oo and 
L — ► oo (with W > L), nor the requirement of a large 



number of propagating modes can be fulfilled (see the 
relevant assumptions in Ref[25]). It is noteworthy that 
the wide end-contacted devices considered here resemble 
experimental setups with invasive contacts, which do not 
show universal behavior, either. 
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